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We present first results from a simulation of quenched overlap fermions with improved gauge field action. Among 
the quantities we study are the spectral properties of the overlap operator, the chiral condensate and topological 
charge, quark and hadron masses, and selected nucleon matrix elements. To make contact with continuum physics, 
we compute the renormalization constants of quark bilinear operators in perturbation theory and beyond. 



1. INTRODUCTION 

Lattice calculations at small quark masses, i.e. 
in the chiral regime, require actions with good 
chiral properties. Overlap fermions PP have an 
exact chiral symmetry on the lattice J2j and thus 
are predestined for this task. A further advantage 
of overlap fermions is that they are automatically 
0(a) improved 

The massive overlap operator is defined by 



D= 1- 



2p 



D 



N 



(1) 



* Talks presented by M. Giirtler, R. Horsley, H. Per It and 
T. Streuer at Lattice 2003. 



D N = P -(l + n ^=),X = D w -P, (2) 

where Dyj is the Wilson-Dirac operator. We as- 
sume r = 1 throughout this paper. The operator 
Djq has n_ + n + exact zero modes, D^ip n = 0, 
n_ (n+) being the number of modes with negative 
(positive) chirality, 75^™ = -ip n {l^n = +ip n )- 
The index of -Djv is thus given by v = ri_ — n + . 
The 'continuous' modes A, Dn4>\ — ^V'a, having 
{iPx,lb4>\) — 0, come in complex conjugate pairs 
A, A*. 

To compute the 'sign function' 



sgn(Y) = 



X 



l5 sgn{H),H = l5 X, (3) 



Vxvc 

we use Zolotarev's optimal rational approxima- 
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tion |IJ. To improve the accuracy of the rational 
approximation, and to reduce the number of iter- 
ations in the inner inversion, we project out the 
■~ 16 lowest eigenvalues of H. The approxima- 
tion of the 'sign function' is done to better than 
5 • 1CT 7 in the interval [0.1, 2.4]. We use a multi- 
mass conjugate gradient solver in both the inner 
and outer inversions. For any given quark mass 
this allows to compute propagators for a whole set 
of higher quark masses at very little extra cost in 
CPU time. For the inner and outer inversions a 
stopping criterion of 10 -6 and 5 • 1CP 6 , respec- 
tively, is employed. 

It is important to use a good gauge field action, 
because the inversion time of the fermion matrix 
is greatly reduced for improved gauge field ac- 
tions. We use the Liischer-Weisz action jS] 



S[U] = \ [c l R ° Tr (! - ^Plaquettc) 

^ plaqucttc 

+ Cl Q ReTr t 1 - Rectangle) (4) 
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£ -ReTr(l-C/ p . 
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with coefficients c\, c-i (cq + 8ci + 8C2 = 1) taken 
from tadpole improved perturbation theory 
In Fig. 1 we compare the condition number of 
the Wilson and tadpole improved Liischer-Weisz 
gauge field action. We see that the condition 
number is a factor of > 3 larger for the Wilson 
action. This is largely due to the fact that the 
Liischer-Weisz action suppresses dislocations [7] 
and thus greatly reduces the number of (unphys- 
ical) zero modes. Topological studies using the 
Wilson gauge field action should be taken with 
caution. 

The calculations are mainly done on the 16 3 32 
lattice at (3 = (6/g 2 )c = 8.45, where ci = 
-0.15486 and c 2 = -0.013407. The correspond- 
ing lattice spacing is a — 0.095 fm if we use 
ro = 0.5 fm to set the scale. We have taken 
p = 1.4, which we have found to be the optimal 
choice. In addition, part of the calculations have 
been done on the 12 3 24 lattice at (3 — 8.1 to test 
for scaling. The lattice spacing here is a — 0.125 



Figure 1. The condition number for the Wilson 
and tadpole improved Liischer-Weisz gauge field 
action on the 16 3 32 lattice at lattice spacing a 
0.1 fm for p = 1.4, as a function of configuration 
number. 

fm [5]. Both lattices have the same physical vol- 
ume V. 

2. SPECTRAL PROPERTIES 

The eigenvalues of Dm lie on a circle around (p,0) 
with radius p in the complex plane. The improved 
operator D™ p = (1 — clDn /2p)~ 1 Dn EH projects 
the eigenvalues of D n stereographically onto the 
imaginary axis. The 'continuous' eigenvalues of 
Z?™ p come in pairs ±iA, while the zero modes 
are untouched. 

We have computed the lowest « 140 eigenvalues 
of Dm on both our lattices for 0(250) gauge field 
configurations each. We employed the Arnoldi 
algorithm as provided by the ARPACK package. 

If one defines the topological charge density by 
q{x) = (l/2)Tr(75Div(a;,a;)) 0, the total charge 
Q is given by the index of Dm'- 



(5) 



Using this definition we have computed the topo- 
logical susceptibility Xtop — (Q 2 )/V and find 



p 


Xtop 


8.10 


(195(5) MeV) 4 


8.45 


(187(5) MeV) 4 



(G) 



The spectral density of the 'continuous' modes is 



3 



formally given by 

p(a) = ^(£*(a-a)>, 



(7) 



where the sum extends over the (nonzero) eigen- 
values ±iA of D™ lp . In our case < |A| < 800 
MeV. In practice one groups the eigenvalues into 
bins, whose size will depend on the statistics. In 
the infinite volume E = —(-ipip) = — 7rp(Q). In the 
finite volume and for small eigenvalues the spec- 
tral density can be computed from the chiral low- 
energy effective theory. For A < Et, Et being the 
Thouless energy Et « f^/Tly/V, the low-energy 
effective partition function is dominated by the 
zero momentum modes, and the zero-momentum 
approximation of the chiral low-energy effective 
theory is equivalent to chiral random matrix the- 
ory. In random matrix theory the microscopic 
spectral density in the sector of fixed topological 
charge Q, 



p^\T,V\) = lim ^p(X) 



1 



(8) 



with YjVX kept finite, is given by |§] 



pf\ x ) = !0Jjo|(*) " J\q\ + i{x)J\q\-i{x)), (9) 



where J n (x) are Bessel functions. Thus, for A < 
Et the spectral density is given by 



P (A) = E^ W (Q)^ Q) (EFA), 



(10) 



where w(Q) is the weight of the sector of topo- 
logical charge Q with J2q w (Q) = 1- Taking 
w(Q) from our 'measured' charge distributions, 
we then may obtain E by fitting (|10fl to our 
data. In Fig. [21 we show the data together with 
the fit for f3 = 8.45. We observe that ran- 
dom matrix theory describes the data well up to 
A sa 150 MeV. For the (unrenormalized) chiral 
condensate we find E = (242(8) MeV) 3 at (3 = 8.1 
and E = (249(9) MeV) 3 at (3 = 8.45. We will give 
renormalized values after we have computed the 
renormalization constants. 

Random matrix theory predicts furthermore the 
distribution of the smallest nonzero eigenvalue in 




Figure 2. The spectral density p(X) at (3 = 8.45, 
together with a fit of random matrix theory to 
the data. Both p and A are given in lattice units. 

the sector of topological charge Q, which we call 
/9^(EVA). The first few expressions are [5]: 



Pse( x ) 
Pse( x ) 
Psk x ) 



-x 2 /2 



(11) 



-x 2 /2 



(ll(x)-h(x)I 3 (x)). 



After having determined E from the spectral den- 
sity of all eigenvalues, these formulae contain no 
free parameters anymore, and thus can be com- 
pared directly to our 'measured' distributions. In 
Fig. 3 we show our data together with the pre- 
dictions of random matrix theory. We find good 
agreement. 

3. RENORMALIZATION 

To obtain physical results from lattice calcula- 
tions of hadron matrix elements the underlying 
operators have to be renormalized. Let us denote 
the lattice regularized operators by 0{a). We 
then define renormalized operators R (p) by in- 
troducing the renormalization constant Zo(ap): 



R {p) = Z (ap)0(a) 



(12) 



The renormalization constant Zo(ap) is found by 
imposing the renormalization condition 



A„ 



= Z i ,(ap)Z (ap)- 1 A t S cc 
+ other Dirac structures, 



(13) 
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Figure 3. Distributions of the smallest nonzero eigenvalue p^(x) as a function of x = T,VX for \Q\ = 
0, • • • , 5 on the 16 3 32 lattice at j3 = 8.45, together with the predictions of random matrix theory, in lattice 
units. 



where Z$ is the wave function renormalization 
constant to be defined below and Aq is the 
forward (quark line) amputated Green function 
computed between off-shell quark states with 4- 
momcntum p. We consider quark bilinear oper- 
ators only. The renormalized operator O r (/j,) is 
independent of the regularization scheme, but will 
depend on the external states and on the gauge. 
The operator matrix elements can be converted 
to more popular schemes like MS by means of 
continuum perturbation theory. 

In the following we present first results of a 
perturbative calculation of renormalization con- 
stants for overlap fermions and improved gauge 
field action. So far results for overlap fermions 
are only known for the Wilson gauge field ac- 
tion [lUlllj . The numbers given below hold for 
ci = -0.15486 and c 2 = -0.013407 (i.e. the im- 
provement coefficients at (3 = 8.45), while we keep 
the (bare) coupling constant g arbitrary, which 
allows us to exchange it for a better expansion 
parameter. 

The lattice Feynman rules for overlap fermions 
have been derived in [12113) . The gluon propa- 
gator for improved gauge field action is known in 



four dimensions only |14I15| . A suitable form for 
dimensional regularization is 

D™*> = D™*°» + AD^, (14) 

where 

is the standard Wilson propagator with = 
(2/a) sm(ak fl /2), and £ specifies the gauge. The 
Landau gauge corresponds to £ = 1, the Feyn- 
man gauge to £ = 0. The six-link interaction 
terms are contained in AD^ = Dffi - D™ Xson . 
While Dj^ llson leads to infrared divergent expres- 
sions, which have to be regularized, expressions 
involving AD^ are infrared finite and thus can 
be computed in four dimensions. We may write 

4 

AZV = ^^D n (fc, Cl ,c 2 )g" 

n=0 

(16) 

m-\-n— 4 
m,n— 
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The coefficient functions D n and D m n — D n m 
are rational functions involving k and ci,2- Ex- 
plicit expressions are given in J^j • Both functions 
vanish in the limit a. 2 — > 0. 

SeZ/ energy 

The inverse of the massless quark propagator can 
be written in the form 



s- 1 -- 

with 



lat 



^lat 



g 2 CF 

16vr 2 



iyEi(oV)- 



For the quark self energy Si we find 
E^aY) = (1-OMaY) 



4.7920U + & S 



with 

6 S = -16.179. 



(17) 
(18) 

(19) 
(20) 



From l|19() we obtain the quark wave function 
renormalization constant 



1 



9 2 C F 
16tt 2 



Ei(aV) 



(21) 



Local operators 

Let us consider local operators 
Ox =rpT x iP = X 
now with 



(22) 





T x 




1 


P 


75 


V 


7m 


A 


7m 75 


T 


CjUi/75 



To find the renormalization constants we have to 
compute the amputated Green functions Ao x = 
A x . We obtain 



A S ' P = {1,75} 



9 2 C F 



16tt 2 

5.7920U+6 S ,p]{1,75}, 



[-(4-0 ln(aV) 



(23) 



= {7„7,75} + f^{h(l-0 In(aV) 



- 4.79201^ + 6y iA ] 7A1 
-2(1-0^} {1,75}, 

A^ = ^75 + ^Kln(aY) 

- 3.7920ie + 6 T ]^ 7 5, 
where 

b s ,p = 10.512, 
6y jA = 6.228, 
6 T = 3.900. 



(24) 



(25) 



(26) 



Using l|13[) we then arrive at the renormalization 
constants 



9 2 G F[ 



Zs.p = 1 ~ t^-tt [-6 ln(o^) - £ + frs.p + &s] , 
I67H 

16tt 



(27) 



1 



S^F 



16tt 2 

In the MS scheme this gives 



[21n(a(i)+£ + b T + bx\. 



9 2 C F 



yMS _ -i 



16tt 



z ~ms 1 5 ^f 



' V ' A ~ * Kirr- 
= 1 



167T 2 

One-link operator 



-6 ln(a/i) - 5 + 6 S: p + & E ] , 
(28) 

2 ln(a/i) + 1 + 6 T + 6 E ] • 



The one-link operator that will be of interest to 
us here is 



O^v = -j^ln Dvip - Traces. 



(29) 



Two different irreducible representations of 0^ v 
under the hypercubic group have been considered 
in the literature |17|: 



O a = C{14}, 

1 3 

Ob = Oa — — ^ On, 



(30) 
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where {• ■ ■} denotes symmetrization of the in- 
dices. Because of space limitations we will not 
give the Green functions here but only state the 
final results. For the renormalization constants of 
the operators O a ,b we obtain in the M S scheme 



7 MS 



with 



= 1- 



9 2 C F 



16tH 



16 



ln(a/z) + 6 a ,6- 



40 



(31) 



b a = -6.516, 
b b = -5.617. 

Improvement 



(32) 



Lattice perturbation theory is known to converge 
badly due to the appearance of (gluon) tadpole 
diagrams, which are lattice artifacts and which 
make the bare coupling g into a poor expansion 
parameter. It was proposed [HQ that the pertur- 
bative series should be rearranged in order to get 
rid of the tadpole contributions. We have made 
use of this observation already in tuning the coef- 
ficients of the Liischer-Weisz action For the 
renormalization constants this rearrangement is 
done in |16| . following a method similar to that 
in Writing 



Z 



o 



1 



(33) 



we arrive at the tadpole improved result 

C F g 2 



7 TI 



P- 



4u 



1 



4 
P 



nnjku j, 



16-7T 2 tig 



Be 



(34) 



Alternatively, for the local operators (with ri£> = 
0) one may improve the perturbative result by 
writing |2*T] 



ryV l ^llOIipClt 



1 



C F g 2 
16tt 2 



B x 



(36) 



where Zy° npcrt is the nonperturbatively deter- 
mined renormalization constant of the local vec- 
tor current and Zy Clt the perturbatively com- 
puted one. A similar procedure can be envis- 
aged for the one- and higher- link operators. In 
that case z v onpert will have to be replaced by the 
appropriate nonperturbatively determined renor- 
malization constant of the one- and higher-link 
operator, respectively. 

In Table 1 we compare the results of the vari- 
ous improvement schemes with the perturbative 
result. For z v ° npeTt we have taken the nonpertur- 
bative result Za = 1.416 derived below (eq. l|4X|l). 






Zq 


7 T1 


Zjq 


S,P 


1.150 


1.190 


1.430 


V,A 


1.140 


1.171 


1.416 


T 


1.159 


1.207 


1.443 


o a 


1.257 


1.335 




o b 


1.244 


1.308 





Table 1 

Comparison of renormalization constants at the 
scale a/i = 1 for /3 
the MS scheme. 



.45 and Z™ npclt = 1.416 in 



where 



UQ — ( gTr t/plaquette) ' 



(35) 



no is the number of covariant derivatives and |20| 
k u = 0.7325 vr 2 . The difference to the Wilson 
gauge field action is that (1 — n^n 2 now has to be 
replaced by (1 — 4/p — n£,)k ui and there is an ad- 
ditional prefactor of p/(p — 4 — 4u ) ■ At our value 
of P (i.e. g 2 = 1.6658) we find u% = 0.65176. 



4. CHIRAL CONDENSATE 

Knowing Z$ , we can now compute the renormal- 
ized chiral condensate 



{^) R (fx) = -Z s (a^)E. 



(37) 



It is traditional to quote numbers for /i = 2 GeV. 
In the MS scheme we obtain Z^(2 GeV) = 1.253 
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Figure 4. Effective pion masses for our four dif- 
ferent quark masses from (At(t)y4.4(0)) . 

at /3 = 8.1 and Zj / (2GeV) = 1.184 at /3 = 8.45. 
This gives 



a 


(ViV) MS (2GeV) 


8.10 


(261(9) MeV) 3 


8.45 


(263(9) MeV) 3 



We find good agreement with scaling. 

At = 8.45 we have Z% z (2GeV) = 1.426. Us- 
ing this value for the renormalization constant we 
obtain (^) JTW (2 GeV) = (280(10) MeV) 3 . 

5. QUARK AND HADRON MASSES 

The results presented in this and the next Section 
refer to the 16 3 32 lattice at j3 = 8.45. The calcu- 
lations are done for four different quark masses, 
am q = 0.028, 0.056, 0.098 and 0.140 (correspond- 
ing to a/x = am q /2p = 0.01, 0.02, 0.035 and 0.05.) 
To increase the overlap of meson and baryon oper- 
ators with the ground state wave function we use 
smeared sources [321 with k s = 0.21 and N s = 50. 
The calculations of two-point functions are based 
on O(100) configurations. The code has partly 
been written in SZIN .22J, which has the advan- 
tage of being flexible and machine independent. 

Pion mass 

We compute the pion mass from various cor- 
relation functions: (P(t)P(0)), (A 4 (t)P(0)} and 



Figure 5. Effective pion masses from 
(^4(^^44(0)) before (solid symbols) and af- 
ter (open symbols) removal of zero modes from 
the quark propagators. 

(^4(^)^4(0)), where P is the pseudoscalar den- 
sity and the axial vector current. In Fig. 4 we 
show the effective mass plot of (Ai(i)At(O)) for 
our four quark masses. This correlation function 
is least affected by finite size corrections induced 
by zero mode contributions 24 25 . The result- 
ing pion masses are given in Table 2. Our lightest 
pion mass is « 400 MeV, which is compatible with 
quenched Wilson fermion calculations. 

The correlation functions (P(t)P(0)) and 
(A4(t)P(0)) give compatible results, albeit with 
larger error bars. If one explicitly removes the 
zero mode contributions from the quark propa- 
gators one arrives at the result shown in Fig. 5. 
It appears that the pion mass obtained from 
(Ai(i)At(O)) increases by w 10% (5%) at the 



am q 




m % [MeV] 


0.028 


0.211(3) 


439(6) 


0.056 


0.288(2) 


599(4) 


0.098 


0.383(2) 


797(4) 


0.140 


0.466(2) 


969(4) 



Table 2 

Pion masses from (Ai(t)A4(0)). We have used 
a = 0.095 fm to convert the lattice numbers to 
physical units. 
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Figure 6. The pion mass (am^) 2 as a function of 
am q , together with the fit. 




c0 

1.4 - 



1.2 

0.05 0.1 0.15 0.2 

am q 

Figure 7. The ratio amj /m q as a function of am q , 
together with the fit. 



lowest (highest) quark mass. Whether this is 
a valid procedure to reduce finite size effects is 
not clear to us, because it involves a nonlocal 
step But it should be taken as a warning 

that finite size effects might still be large. 

In Fig. 6 we plot m 2 against m qi and in Fig. 7 
we show the deviation of m 2 from linearity. 
Quenched chiral perturbation theory predicts, in 
the infinite volume, 

ml = Am q (l - 5[\n(Am q /A 2 x ) + 1]) + 0(m 2 ). (39) 
We fit our data by 

ml = Am q + Bm q \nm q + Cm 2 q . (40) 

The result of the fit is shown by the curves in 
Figs. 6 and 7. Using A x = 4irU (U = 93 MeV), 
we derive S — 0.26(9). This number agrees, 
within error bars, with what one would expect. 
However, before one can draw any conlusions, a 
careful study of finite size effects must be per- 
formed. 

Quark masses 

We determine the bare light and strange quark 
masses, mi — (m u + md)/2 and m s , from the 



physical pion and kaon masses |27I28| : 

m T P h ys = Ami + B m£\nm e + C mj, (41) 

m ir P hys = A(mi + m s )/2 + B (mi + m s )/2 
x / m s lnm s - m e Inmi _ \ 
V m s — mi ) 
+ C[{mi + m s )/2} 2 . 

We find ami = 0.0020(3) and am s = 0.068(2). 
The renormalized quark mass is given by 
= Z m (afi)m q with Z m = l/Z s . Using 
Zj J (2GeV) = 1.184 (and a = 0.095 fm), we ob- 
tain 

mf^{2 GeV) = 3.5(3) MeV, (43) 
raf^(2 GeV) = 119(4) MeV. (44) 

If we use Z|' / (2GeV) = 1.426 instead, the num- 
bers reduce to 

mf^(2 GeV) = 2.9(3) MeV, (45) 
mf^(2GeV) = 99(3) MeV. (46) 

The light quark mass mi should not be taken too 
seriously, because it is strongly affected by chiral 
logarithms. 
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Figure 8. The renormalization constant Za as a 
function of am q , together with the chiral extrap- 
olation. 



Figure 9. The nucleon mass (toajto) 2 as function 
of am q , together with the experimental result (*), 
using ro/a = 5.26. 



Nonperturbative determination of Za 



Nucleon mass 



The renormalization constant Za can be com- 
puted nonperturbatively from the Ward identity 

(cf. m 



t^oo mv {A±{t)P{Q)) 



(47) 



The result is plotted in Fig. 8. We see that Za 
depends only weakly on the quark mass. A linear 
extrapolation to the chiral limit gives 



Z A = 1.416(2). 



(48) 



This number lies 20% above the tadpole improved 
perturbative result. 



am q 




m N [McV] 


0.028 


0.60(3) 


1250(60) 


0.056 


0.65(1) 


1350(20) 


0.098 


0.775(10) 


1612(20) 


0.140 


0.875(5) 


1820(10) 



Table 3 

Nucleon masses. We have used a — 0.095 fm to 
convert the lattice numbers to physical units. 



Our results for the nucleon mass are given in Ta- 
ble 3. In Fig. 9 we plot {tun To) 2 against ra q . 
Except for the lowest mass, which currently is 
not very well determined, the data points lie on a 
straight line, which extrapolates surprisingly well 
to the experimental value. 

6. NUCLEON MATRIX ELEMENTS 

While hadron masses are determined from two- 
point correlation functions, nucleon matrix ele- 
ments of quark bilinear operators 0, 



(N\0\N), (N\N)=2m N , 



(49) 



are derived from ratios of three-point to two-point 
functions, 



R = 



(N(t)O(r)N(0)) 1 
(N(t)N(0)) 2m N 



(N\0\N), (50) 



which for t 3> r ^> are proportional to the de- 
sired matrix element, as shown on the r.h.s. of 
eq. l|5Ujl. Here N(t) is a suitable baryon operator. 
We consider nucleons of zero momentum only. 
This technique is by now standard. For details 
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Figure 10. The ratio R plotted against the posi- 
tion t of the operator Ob for am q = 0.14. The 
source (N) is placed at time slice t = 0, and the 
sink (N) at t = 13. From a plateau at t 3> t ^> 
we can determine R. 

and the numerical implementation the reader is 
referred to |22I29| . The calculation of three-point 
functions requires additional inversions and there- 
fore is computationally expensive. The calcula- 
tions in this Section are based on O(50) configu- 
rations, so the results must be regarded as very 
preliminary. 

While hadron masses are automatically 0(a) im- 
proved, operator matrix elements are generally 
not, except in the chiral limit. However, in dis- 
tinction to improved Wilson fermions, improve- 
ment for overlap fermions is very simple and uni- 
versal 0: 

o imp = o. (51) 

This a major advantage, as we do not have to 
compute matrix elements of higher dimensional 
operators, nor their associated improvement co- 
efficients. 

We shall consider two operators, the local vector 
current and the one-link operator Ob- The 
nuclcon matrix element of Ob yields the first mo- 
ment of the unpolarized nucleon structure func- 
tion. The results presented in this Section refer 
to improved, nonsinglet operators. 



A typical ratio R is shown in Fig. 10 for the oper- 
ator Ob from which we can find the bare nucleon 
matrix element. 

Local vector current 

Let us first look at the local vector current. Due 
to charge conservation 

(N\Vf\N) = Z v (NIV^N) = 2m N , (52) 

which may be used to compute Zy |30| . We ex- 
pect Zy = Za, and the idea is to test this hypoth- 
esis. In Fig. 11 we plot Z v against am q . We com- 
pare this result with Za computed from eq. (|47(l 
in Section 5. We see that Zy approaches Za in 
the chiral limit, while at our largest quark mass 
the numbers differ by a few percent. Though the 
operators are 0(a) improved, discretization errors 
~ («m T ) 2 are still possible, which might explain 
the small slope. A linear extrapolation to the chi- 
ral limit gives 

Z v = 1.426(7). (53) 

First moment of the structure junction 

Let us now turn to the operator Ob- The nucleon 
matrix element of the operator O a is much harder 
to compute, as it requires a nonzero nucleon mo- 
mentum. The matrix element of the renormalized 
operator Ob gives [221 

(N\0?\N) = Z b (N\O b \N) 

(54) 

= -2m 2 N (x). 

For the tadpole improved renormalization con- 
stant in the ~MS scheme we find Z h MS (2 GcV) = 
1.314. In Fig. 12 we show the first moment of 
the nonsinglet nucleon structure function (x) in 
the MS scheme at fi = 2 GeV. The numbers are 
in surprisingly good agreement with the experi- 
mental value. No nonanalytic behavior is needed 
to reconcile experimental and theoretical results. 
Previous calculations using Wilson fermions gave 
much larger values. 
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Figure 11. The rcnormalization constant Zy from 
this calculation (solid symbols) compared with 
Za from Section 5 (open symbols) as a function 
of am q , together with linear extrapolations to the 
chiral limit. The symbols at the lowest quark 
mass are displaced by a small amount so that they 
do not overlap. 



7. SUMMARY 

Overlap fcrmions have many advantages over 
Wilson and staggered fermions. They provide 
an implementation of lattice fermions with ex- 
act chiral symmetry, even at finite lattice spac- 
ing. In addition, they are automatically 0(a) im- 
proved, and the task of operator renormalization 
is greatly reduced. 

Calculations with overlap fermions on fine 
grained, phenomenologically relevant lattices are 
progressing rapidly. They are computationally 
costly, but by using an improved gauge field ac- 
tion and projecting out the lowest lying eigenval- 
ues the condition number can be substantially re- 
duced. We have tested the predictions of random 
matrix theory and find good agreement with the 
unfolded distributions of the smallest eigenvalue 
in topological sectors up to \Q\ = 5. The renor- 
malization constants of local and one-link oper- 
ators have been computed perturbatively for the 
tadpole improved Liischer-Weisz gauge field ac- 
tion. The one-loop corrections turn out to be rela- 



Figure 12. The first moment of the nucleon struc- 
ture function (x) at /i = 2 GeV as a function of 
am q . Also shown is the phenomenological value 

(*)• 

tively large, and in most cases significantly larger 
than in the improved Wilson fermion case |19| . 
which calls for a nonperturbative determination. 
As a first step, we have computed Za and Zy non- 
perturbatively. We find that Za and Zy agree 
in the chiral limit, as expected. The nonper- 
turbative numbers turn out to lie 20% higher 
than the tadpole improved perturbative values 
(at a = 0.095 fm). On the phenomenological side 
we have computed the topological susceptibility 
and the chiral condensate from the spectrum of 
low-lying eigenvalues. The topological suscepti- 
bility is found to be in good agreement with the 
Witten-Veneziano formula 

m v' = —ffxtop, (55) 

giving m v ' — 920(50) MeV. We employed random 
matrix theory to derive the chiral condensate in 
the infinite volume. We went on to compute the 
pion, nucleon and quark masses. We find some 
signal for chiral logarithms. Both r$/a 6 , using 
r = 0.5 fm, and the nucleon mass give compati- 
ble values for the lattice spacing. Our lowest pion 
mass so far is « 400 MeV. There is quite a big 
uncertainty in the strange quark mass due to an 
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uncertainty in Z m . Finally, we were able to com- 
pute the first moment of the nucleon structure 
function. Good agreement with the experimental 
value has been found. 
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